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Abstract 

We study the partial differential equation 

max{Lu - f,H(Du)} = 

where u is the unknown function, L is a second-order elliptic operator, / is a given 
smooth function and H is a convex function. This is a model equation for Hamilton- 
Jacobi-Bellman equations arising in stochastic singular control. We establish the exis- 
tence of a unique viscosity solution of the Dirichlet problem that has a Holder continu- 
ous gradient. We also show that if H is uniformly convex, the gradient of this solution 
is Lipschitz continuous. 
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1 Introduction 



In this paper, we consider PDE associated with a general class of stochastic singular control 
problems. This is a class of nonlinear, second-order PDE that each have a free boundary 
determined by a convex gradient constraint. Using PDE methods, we show that the Dirichlet 
problem has a unique solution and derive some regularity properties of the gradient of this 
solution. Namely, we establish that the gradient is Holder continuous and also that if H is 
uniformly convex, the gradient is Lipschitz continuous. Finally, we give a brief discussion of 
how this type of equation arises in control theory and show how our regularity results apply 
to the motivating control problems. 
The PDE we focus on is 



where O C M n is open and bounded with smooth boundary 80 and / is a smooth, non- 
negative function on O. We assume that L is the linear differential operator 

Li/j(x) := -a(x) ■ D 2 tfj + b(x) ■ Dtf; + c(x)^, V e C 2 (0) 

with smooth coefficients a : O S(n), b : O — > lR n and c : O — > R. Here A - B := tr A f B and 
S(n) denotes the set of n x n symmetric matrices with real entries. We shall further assume 
that L is (uniformly) elliptic: 



max{Z/u - /, H(Du)} = 0, x G O 

u = 0, x e 80 ' 



(1.1) 



a(x)£-£ > 7|£| 2 , for all xeO^eW 
for some 7 > 0. The final assumption on L that we will make is 



n 



(1.2) 



c(x) >5, x e O 



where 6 is a positive constant. 
Our central result is 



Theorem 1.1. (2) Assume that H : W 1 — > R satisfies 




(1.3) 



Then there is a unique viscosity solution 



«6 4(0)n^(0) 



of (HU), for anyaE (0,1). 



(ii) If, in addition to f 1 1 . 3 [) . H satisfies 



H is uniformly convex 

D 2 H e 



(1.4) 
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thenueC^O). 



We employ techniques from the theory of viscosity solutions of scalar non-linear elliptic 
PDE to prove the existence and uniqueness of solutions of (11.11) . We also use a penalization 
technique similar to the one introduced by L.C. Evans in jl] and refined by M. Wiegner [12] 
and H. Ishii et. al. [8] to establish regularity of solutions. More precisely, we study the 
penalized equation 

Lu* + /3 e (H(Du e )) = f, xeO 

u 6 = 0, x G dO ' 1 ' ' 

where (/3 e ) e >o is what we call a penalty function. /3 e can be thought of as a smoothing of 
z H- (z/e) + . We show that this equation has a unique solution u e that is bounded in W^(0) 
for all e positive and small enough and that u e converges to u in Cl oc (0), though a sequence 
of e tending to 0. 

Our result is novel in the fact that it only assumes convexity of the gradient constraint 
function. Some previous regularity results apply to the case of the gradient constraint 
function 

H(p,x) = \p\-g(x), (p 1 x)eM n x0 1 (1.6) 

where g is a smooth, positive function on O. While we do not consider gradient constraints 
that depend on x, a close inspection of the methods we employ indicate they would apply to 
a large class of these gradient constraints. For instance, if we assume that H G C(M n x O) 
satisfies 

p i— )■ H(p, x) is convex for each x G O 
H(0, x) < for each x G O , 
x M- H(p,x) G C ' 1 ^) 

then a modification of our methods would establish part (i) of Theorem 11.11 for solutions of 
equation (jl. ip with this H. 

Our goal was to identify general structural conditions on the type of gradient constraints 
for which penalization methods are successful at yielding regularity results. We believe that 
assumptions (11. 3p and (11. 4p accomplish this. In addition, we remark that our main regularity 
result involves uniform convexity of H while previous results did not explicitly make this 
assumption. This is because uniform convexity was built in by their choices of gradient 
constraints. For example, we can replace the gradient constraint given by H defined in (11. 6p 
with \Du\ 2 — g(x) 2 < 0, which is uniformly convex in Du. 

It also should be noted that C 2 regularity of solutions has been obtained in the case 
O = M. n and H(jp) — \p\ — 1 by assuming n = 1 or 2, very special structural conditions on 
/ and/or that L has constant coefficients [TUj [TTj . The purpose of this work was to study 
solutions of ( 11. ip on a general bounded domain in MJ 1 and estimate solutions for a general 
gradient constraint function and a general elliptic operator. 

The organization of this paper is as follows. In section [21 we show that ( 11. ip has a unique 
viscosity solution by establishing a comparison principle for sub and supersolutions. Then we 
pursue the regularity of this solution in section [3] by studying the penalized equation ( 11. 5p . 
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In the final section, we use our uniform estimates to pass to the limit as e — > + and prove 
Theorem 11.11 Before performing our analysis, let us discuss the motivating applications in 
singular control theory and give a probabilistic interpretation of solutions of (11. ip . 

Probabilistic interpretation of solutions. Assume that (Q, J 7 , P) is a probability 
space equipped with a standard n-dimensional Brownian motion (W(t),t > 0). A control 
process is a pair (p, £) such that 

'(p(i),£(*))eR n xR, t >0, 
(p, £) is adapted to the filtration generated by 
\p(t)\ = l,t> 0,a.s. 

,£(0) = 0, t y is non- decreasing and is left continuous with right hand limits a.s. 
Now, let I be the support function of a nonempty, closed, convex set K C lR n . That is, 

£(v) := supt; -p, veR n . (1.7) 

We consider the stochastic control problem 

u{x) : = infE x f e -fi <*"•*(.'))*> [f(X p ^(t))dt + £(p(t))d£(t)] , z G O. (1.8) 

Here X p ^ satisfies the stochastic differential equation (SDE) 

dX(t) = -b(X(t))dt + a(X(t))dW(t) - p(t)d£(t), t > 
X(0) = x 

and r = inf{t > : X p ^[t) £ O}. We are assuming that o~,b,c are smooth on O and that 
the above SDE has as unique solution (in law) for each x G O and control process (p, £). 
In general, X will not have continuous sample paths and so it is regarded as a "singularly" 
controlled process. Therefore, we say that u is the value function of a problem of stochastic 
singular control. 

W. Fleming and H. Soner have shown that if the value function u satisfies a natural 
dynamic programming principle, then u is a viscosity solution of a Hamilton- Jacobi-Bellman 
(HJB) equation of the form (II. ip (Theorem 5.1, section VIII. 5 in [6]). This result, restated 
below, provides the connection between equation (11. ip and stochastic singular control. 

Theorem 1.2. Assume that for each stopping time 9 (with respect to the filtration generated 
by W) and x G O, 

u(x) = inf E x \e~ ^ W^^X^t A 6)) 

+ / e~ J"" c ( XP ' 5 ( s )) ds [/(X^(t))dt + £(p(t))^(t)] 
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Then the value function u is a viscosity solution of HJB equation 

max {-\a{x)a\x) ■ D 2 u + b(x) ■ Du + c(x)u - f(x),H(Du)} = 0, x E O 



where 



(1.9) 
u = 0, x G dO 



H(p) = max {p- v-£(v)}, pel". (1.10) 



In particular, u is a viscosity solution of (11.11) with 



a(x) := -a(x)a t (x), x e O. (1.11} 



In view of Theorem ll.il we have the following corollary which asserts that singular stochastic 
control problems as described above have C 1 value functions. 

Corollary 1.3. Assume the hypotheses of Theorem that a : O — >■ S(n) defined in (11. lip 
satisfies the uniform ellipticity condition (jl.2p , and that G M n an interior point of the 
convex set K associated with i. 

(i) Then u given by (11.81) is the unique viscosity solution of (II. 9p . In particular, u G C]^{0) 
for each < a < 1. 

(ii) Further suppose i is the support function of K := {p G K n : G(p) < 0} where G satisfies 
f Oj) and (HU). T/ien m G (7^(0). 



Proof, (i) If lies in the interior of K, B$(0) C K for some 5 > 0. It is immediate from 
(ll.7p that £ > 5 and in particular that H(0) < —6 < 0. The conclusion follows at once from 
Theorem 11.11 

(ii) A standard fact about the support function £ of a closed, convex subset K C M n is 
that 

p| { P 6l n :t)-p< £(u)}. 

(Theorem 8.24 in [9]). Therefore, we have by hypothesis 

{pel": G(p) < 0} = p) {p G M n : v ■ p < £(v)} . 

\v\=l 

In view of (ll.lOp . it is also plain that 

{p G W n : H{p) < 0} = P {p6R":t)-p< £(w)} . 



v=l 



Consequently, u solves the PDE 
1 



max 



2 a{x)a t {x) ■ D 2 u + b(x) ■ Du + c(x)u - f{x), G(Du) } = 0, x G O. 



From the assumptions made on G, it follows from Theorem 1 1 . 1 1 1 hat u G C lo ' c (0). □ 
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2 Comparison principle 



In this section, we will verify a fundamental comparison principle among viscosity sub- and 
supersolutions of the PDE 

max{Lu - f,H(Du)} = 0. 

Although it is now well known how to establish comparison principles of elliptic nonlinear 
PDE, the problem we have must be individually addressed because of the gradient constraint. 
In what follows, we use the notation and several basic results found in standard sources on 
viscosity solutions such as [U [2] and the references therein. Throughout this paper, all PDE 
and partial differential inequalities will be interpreted in the viscosity sense. Therefore, we 
may sometimes omit the term "viscosity" when we mention solutions, subsolutions, and 
supersolutions. 

Proposition 2.1. Assume u is a subsolution of ( II. ip and v is a supersolution of (11. ip . If 

u<v ondO and u G L°°(dO), 

then u < v in O. 

Formal Proof. Before proving the above proposition, we give a formal proof (i.e. assuming 
u, v G C 2 (0)) that will help motivate a rigorous argument. Fix e G (0, 1) and set 

w e (x) = eu(x) — v(x), x G O. 
w e G USC(O) and thus achieves its maximum at some x t G O. If x t G dO, then 

w e (x e ) = -(1 - e)u(x e ) + u(x e ) - v(x e ) < -(1 - e)u(x e ) < (1 - e)|«U°°(0o)- 
If x e G O, then by calculus 

= Dw t (x e ) = eDu(x e ) — Dv(x t ) 
> D 2 w e (x t ) = eD 2 u(x t ) - D 2 v(x e ) 

As H(Du{x e )) < 0, 

H(Dv(x e )) = H(eDu(x e )) < eH(Du(x e )) + (1 - e)ff(0) < (2.1) 

by (II. 3p . In particular, since v is a supersolution, we have that 

Lv{x e ) - f(x e ) > 0. 

Therefore, 

c(x e )w(x t ) < L(eu — v)(x e ) 

< -(l- e )f(x e ) 

< 

and hence w e (x e ) < 0. In either case, w e < C(l — e), and letting e — > 1 _ gives u < v. 
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□ 

Proof, (of the proposition) 1. Fix e G (0, 1) and set 

w' n {x, y) = eu(x) - v(y) - — \x - y\ 2 , x,y eO 

for r) > 0. w r ' G U SC(0 x O) and so has a maximum at some point (x^, y n ) 6 0xO. As O 
is compact, (x v ,y v ) has a limit point of the form (x e ,x e ) through some sequence of 77 — > + , 
where x t is a maximizing point of x 1— > eu(x) — v(x), and through this sequence of 77 — > + 



77 



2 

- (2.2) 



(Lemma 3.1 in [2]). 

2. If a; e G 90, we have from the definition of w v and our assumptions that 

eu(x) — v(x) < eu(x e ) — v(x e ) = —(1 — e)u(x e ) + u(x e ) — v(x e ) < C(l — e), x G O. 

Now we assume that x e G O and without any loss of generality that (x v ,y n ) G O x O for 
77 > 0. According to the Theorem of Sums (Theorem 3.2 in [2]), for each p > there are 
X, Y G S (n) such that 

V 



V 



XjeJ (eu)(x v 
y) G J 2 '-v(y v ) 



and 



Here 



X 
-F 



A = D 



2\ x -y\ 



2'/ 



<A + pA\ (2.3) 



1/7 -7 
-7 7 



x=x v ,y=y v T] 

and 7 is the n x n identity matrix. In particular, choosing p = 77 in ( 12.31) implies the matrix 
inequality 

-°,H(v/)- 

3. Since 77 is a viscosity subsolution 

max j-a(z„) • * + &(*„) • + c(zXz„) - 77 (^pj } < 0, (2.5) 

and since 77 is a viscosity supersolution 

max |-o(y„) • F + • + c(i/„Mi/„) - f(y v ), 77 } > 0. (2.6) 



As H {^^j < 0, we have by (OJ) that 

H (^) = H (f*^) < eH (^) + (1 - e)H(0) < 0. 



By (ESI), 

a(y v ) ■ Y + b(x v ) ■ + c(y v )v(y ri ) ~ f(y«) > 0. (2.7) 

Combining (12.5ft and ( 12. 7ft gives, 

ec(x„)u(x^) - c(y v )v{y v ) < a(x n ) ■ X - a{y v ) ■ Y + (b(x v ) - b(y v )) ■ -5 2 



+ef{x v ) - /(l/r?) 



.1 , 



Vrt\ 



< a(x n ) ■ X - a(y v ) ■ Y + Lip (6) — - tfr " + Lip(/)|x^ - 

(2.8) 

Note x H- a 1//2 (x) (the unique positive square root of a(x)) is Lipschitz continuous since 
x \-> a(x) is Lipschitz and a > 7 > 0; indeed 

Lip^) < 

Also note that the 2n x 2n matrix 



is non-negative definite, and by f)2.4p 



a(x v ) ■ X - a(y v ) ■ Y 



By flm 



< 



< 



tr 



< tr 



tr [a(x v )X - a(y v )Y] 

a^ix^ix,) a 1 /2( Xr? ) a V2(^) 

a^{y v )a^{x v ) a^{y,)a^{y v ) 

»WW a 1 /2(^) a V2( ?/j? ) 

a 1 /2( 2/j))a l/2 (X}?) a l/2 ( ^ )a l/2 (2/r?) 



-tr f(a 1/2 ( 
?7 

3Lip(a) 2 



-a 1 /2( 2/r; ))(( a i/2 (Xr;) _ a i/2 ( ^ )) ] 



2Y 



2'/ 



ec 1 



Im(^) - c{y v )v{y v ) < ( ^0^ 



2Lip(6)J 



X 
-Y 

3 / 
7/ 



2// 



Lip(/)|^-^|. (2.9) 
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4. Let (x e ,x e ) be a limit point of {x^^y^) through as sequence of 77 — > + . If x t G dO, we 
have from our remarks above that 

eu(x e ) — v(x e ) < C(l — e). 

If x e G O, we let 77 — > + through the appropriate subsequence in (I2.9P and recall ( 12. 2 p to 
arrive at 

c(x e )(eu(x t ) — v(x e )) < 0. 
This inequality implies eu(x e ) — v(x e ) < 0, and so in either case, 

eu(x) — v(x) < eu(x e ) — v(x e ) < C(l — e), x G O. 

We conclude by letting e — > 1~. □ 

Remark 2.2. With a comparison principle in hand, we can now employ a routine application 
of Perron's method to obtain the existence of solutions. Indeed, observe that 

U EE 



(2.10) 



is a subsolution of (11. IP ; and w, the unique solution of 

Lv = f, x G O 
v = 0, x G 90 

is a supersolution of (11.11) . Therefore, Perron's method applies from which we conclude 

u(x) := sup {w(x) : u< w <u, w is a subsolution of ( II. ip } 

is a viscosity solution (Theorem 4.1 in [2]). By the comparison principle, this solution is 
unique. For the remainder of this paper, we pursue the regularity of solutions of (II. ip . 



3 Penalization method 

In this section, we analyze solutions of the penalized equation (11.51) 

' Lu e + f3 e (H(Du e )) = f, xeO 

u e = 0, xedO' 

where (/3 e ) e >o is a family of functions (/3 e ) e >o satisfying 

'& G C^iR) 

Pe = 0, Z < 

(3 e > 0, z > 
#>0 

[/3 e (z) = z>2e 



(3.1) 
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For each e > 0, we think of e as a type of smoothing of z i— > (z/e) + ; for small e, we think 
of (3 e as a smooth approximation of the set valued mapping 



A>(*) 




t < o 

t = 



It will be a standing assumption that such a family of functions satisfying (13. ip exists. 
The reason for using this approximation is from the following intuition. Since the values 
of f3 E (H(Du e )) can be large when H(Du e ) > and e small, solutions will seek to satisfy 
H(Du e ) < and, in this sense, become closer to satisfying equation (11. ip . 
In our analysis of (jl.5p . we make the following special assumptions on H 

{H G C 2 (R n ) 
D 2 H > 9, for some 9 G (0, 1) . (3.2) 
sup p6R „ \D 2 H(p)\ < oo. 

Admittedly, Theorem 11.11 addresses a much larger class of gradient constraints H. However, 
we will see in the next section that the assumptions made above can be relaxed by replacing 
a general H with (a smoothing of) an appropriate inf-convolution. 

Notice that (I1.5P is a semi-linear, uniformly elliptic PDE with smooth coefficients. By 
our growth assumptions on /3 e and (13. 2p f3 t (H(p)) grows at most quadratically as \p\ — > oo, 
for each e > 0. It follows that (II. 5p has a unique, classical solution u e (Theorem 15.10 in [7]). 
Our goal is to derive W 2,p estimates (1 < p < oo) on solutions that are independent of all 
e > and small. Such estimates would aid us in proving that a subsequence of u e converges 
to u, the solution of (II. ip . in Cl oc {0) as e — > + . To this end, our main result concerning 
the penalization method is as follows. 

Proposition 3.1. Assume that O' GG O, 1 < p < oo and < e < 9. We have the following 
estimates: 

\u e \ w/i,°°(o) ^ C 

for a universal constant C ; 
(ii) 

\u £ \w 2 -p(o') < Ci, 

for some C\ depending on p and O' ; 
(in) 

for some C2 depending on 9, O' , and 

max \D 2 H(p)\. 

IpI<c' 

The proof Proposition 13 . 1 1 is accomplished through the following sequence of lemmas, and 
we will obtain the desired estimates by employing the Bernstein method. 



10 



Lemma 3.2. There is a constant C such that 

\u e (x)\<C, xeO 

for e > 0. 

Proof. Let u be the unique smooth solution of (I2.10p . As u is a supersolution of equation 
(11.51). w e < u\ while u e > 0, since u : x H- is a subsolution of (11.51). Hence, < w e < u. □ 



An immediate corollary of the above proof is 
Corollary 3.3. There is a constant C such that 

\Du e {x)\ < C, xedO 

for e > 0. 

Proof. By the proof of the previous lemma, we have < u e < u with equality holding on 
30. Thus, 

?!<^!<0, xedO 
ov ov 

where v is the outward normal on 80 (which is assumed to be smooth). □ 
Lemma 3.4. There is a constant C such that 

\Du e (x)\<C, xeO 

for0<e<9. 

Proof. 1. It suffices to bound the function 

v € (x) = \Du e (x)\ 2 - \u e (x), XEO 

from above, for some universal (that is, e independent) constant A > 0. To this end, we 
suppress e dependence, function arguments and make use of equation ( II. 5ft to compute the 
following identity 

n 

a ■ D 2 v = 2aD 2 u ■ D 2 u - 2 u XiX] Du ■ Da^ + Dv ■ {(3'DH + b) 

+ 2Du ■ (uDc + cDu + Db-Du- Df) - \{cu + b ■ Du - f) 

+ \(P'Du- DH - P). (3.3) 

2. Since H is uniformly convex with D 2 H > 9 and H(0) < 0, 

p ■ DH(p) - H(p) > e\p\ 2 /2, p G R n . 
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Also note that, /3 e (z) < z(3' e (z) for all z G M which implies 

/3' Du -DH-/3> f3'(Du ■ DH - H) > f3' -\Du\ 2 

Using this inequality, the previous lemma and the uniform ellipticity of a fll.2jl we have from 
( 13. 3 p the estimate 

a-D 2 v> -C\Du\ 2 -C + Dv ((3'DH + b) + ^\Du\ 2 . (3.4) 

3. Let x G O be a maximizing point for v. If xq G dO, a bound on \Du(xq)\ 2 that is 
independent of e > is immediate from the previous corollary. If x G O, then 

Dv{x ) = 0, a(x ) ■ D 2 v{x ) < 0. 

Now, if 0' = (3'(H(Du(x ))) < 1/6 < 1/e, then (3 = (3(H(Du(x ))) < 1 by j3HJ). In 
particular, H(Du{x$)) < 2e < 26 < 2 which implies a bound on |Dw(xo)| independent of 
e G (0,0). If f3'(H(Du(x ))) > 1/6, ^ gives 

0>-C| J Dn(a: )| 2 -C , + ^| J Dn(xo)| 2 , 

which implies a bound on \Du(xo)\ 2 independent of e G (0,6), for A > chosen large 
enough. □ 

Lemma 3.5. Let O' CC O . There is a constant C\ defending on O' such that 

< /3 e (H(Du e (x))) <d, x G Of 

for < e < 6. 

Remark 3.6. To simplify the arguments given below, we assume 

6 = 0, and c = 5 > 0. 

It is straightforward to verify that incorporating more general coefficients b and c is merely 
technical and no new issues arise. 

Proof. 1. It suffices to bound 

v e (x) := r](x)f3 e (H(Du e (x))) + -\Du e (x)\ 2 , x G O 

for each rj G C^°(0), < rj < 1. As before, we will omit the e dependence of u e and v € , omit 
the arguments of functions, write (3 for (3 € (H(Du e )) and compute the following identity for 
a-D 2 v 
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a ■ D 2 v = (a ■ D 2 r])r] + 0'DH ■ Dv + /3"aD 2 uDH ■ D 2 uDH 

/?' j?7 i^a ■ D 2 uD 2 HD 2 u + DH ■ D(5u - f) - ^ u XlXj D aij ■ DH 



} 



- (3DH ■ Dr) - DH ■ D 2 uDu + D l uDu ■ Du 

n 

+ aD 2 u ■ D 2 u + Du ■ D{5u - f) - J] u XiXj D aij ■ Du. (3.5) 

Below, we will make use of the following inequalities 

( /3"aD 2 uDH ■ D 2 uDH > 



< 



a ■ D 2 uD 2 HD 2 u > ^9\D 2 u\ 2 
aD 2 u ■ D 2 u > -f\D 2 u\ 2 
[0 < (3 < C{1 + \D 2 u\} 



(3.6) 



which follow from PDE (11.51) and our various assumptions on /3, a and H (13. ip . (I1.2p . (13. 2p . 

2. Let xo G O be a maximizing point for v. We may as well assume that xq G O; 
otherwise, v < \Du(xq\ 2 /2 < C. In this case, we have 

Dv{xq) = and > a(x ) ■ D 2 v(x Q ). 

From these inequalities, (13. 5 p and (13. 6p . it is straightforward to derive the following inequality 

0> ^\D 2 u\ 2 -C -Cf3' {\D 2 u\ + l} (3.7) 

where C denotes various constants that are independent of e G (0, 9) (but may depend on 
t]{xq), Dt](xq), D 2 7](xq)). All functions in (13 .7p and in the rest of this proof are evaluated at 
x . 

If 

P'<l<l/e or ^\D 2 u\ 2 - C < 0, 

then we have a uniform upper bound on = /3(H(Du(xq))) and also the desired upper 
bound on v (xq). Otherwise, (13. 7p implies 

> i 7 |D 2 u | 2 -C-C(|£> 2 ti| + l) 

and in particular \D 2 u(xo)\ < C. By the last inequality in (I3.6p . v < C\ for some C\ 
independent of e G (0, 6) and only depending on O' . □ 

Thus far, we have established part (i) of Propostion l3.1j the conclusion follows at once from 
the previous lemma. The previous lemma will also aid us in obtaining a pointwise bound on 
the second derivatives of u E and thereby establish part (ii) of Propostion 13.11 To this end, 
we adapt the approach by M. Wiegner [T2] . 
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Lemma 3.7. Let O' CC O and C > as in Lemma 3.4 ■ There is a constant C2 depending 
on 9, O' , max| p |<c \D 2 H{p)\ such that 

|w e |vi/2,cx)( /) < C*2 

for each < e < 9. 

Proof. 1. It is sufficient to bound, for each 77 e C£°(0) with < 77 < 1, the quantity 

M e := max{r](x)\D 2 u e (x)\} 

for all < e < 9. With this in mind, we shall bound the related function 



V '( x ) = -r)(x) 2 \D 2 u e (x)\ 2 + \r](x)P(H(Du e (x))) + ^\Du e (x)\ 2 , xeO 

from above. Here A, fi are constants that will be chosen below. As in previous proofs, we 
shall omit the e dependence of u e ,v e and their derivatives and many times we will write (3 
for /3 e (H(Du £ )). 

Direct computation gives 

r]r] Xi \D 2 u\ 2 + r] 2 D 2 u ■ D 2 u Xi + \(r) Xi /3 + rj^'DH ■ Du Xi ) + fiDu ■ Du x% 
= (vVxixj + r) Xi r) Xj )\D 2 u\ 2 + 2r]r] Xi D 2 u ■ D 2 u Xj + 2f]n X:j D 2 u ■ D 2 u Xi + 
rf{D 2 u Xi ■ D 2 u Xj + D 2 u ■ D 2 u XiXj ) + X [r] XiX fi + r] x J'DH ■ Du Xj + VDH ■ Du Xi 
n { ((3"(DH ■ Du Xi )(DH ■ Du Xj ) + (3'(D 2 HDu Xi ■ Du Xj + DH ■ Du XiX .)) }] 
+/i {Du Xi ■ Du Xj + Du ■ Du XiX . } 

Using the above expressions for v Xi and v XiXj and the fact that u solves the PDE (11. 5ft . it is 
straightforward to verify the following identity 

n 

a ■ D 2 v = (77a • D 2 r] + aDr] ■ r])\D 2 u\ 2 + 4 ^ aijr] Xi D 2 u ■ r]D 2 u X] + 

1,3=1 

n 

aD\ ■ D\ + D 2 u ■ D 2 (6u - /) - ^ u XhXt [2(a Xk ■ D 2 u Xl ) + (a XkXl ■ D 2 u 

v k,l=l 

+li |glD 2 u ■ D 2 u + Du ■ D(5u - f) - u XiXj Daij ■ Du\ + A(a ■ D 2 r])f3 
+/3" (r] 2 D 2 u{D 2 uDH) ■ (D 2 uDH) + ri(3"\-f\D 2 uDH\ 2 ) + f3'DH ■ Dv + 



0' 



A ( aDr] ■ D 2 uDH + 77 <j a ■ D 2 uD 2 HD 2 u + DH ■ D(5u - f) - ^ u XiXj D aij ■ DH 



(3Drj -DH) - r]\D 2 u\ 2 DH ■ Dr] + tfD l u ■ D z uD 2 HD 2 u 



2 n2„ 



(3 
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2. Let Xq be a maximizing point for v. If x G dO, then t> < v(x ) < fi\Du(x )\ 2 < C. 
This of course implies M 2 < v(x ) < C as desired. Now suppose that x G O, so that 

Dv{xq) = and a(xo) ■ D 2 v{xq) < 0. 

Evaluating identity (13. 8 p at xo, employing the assumed matrix inequalities a > 7 (11. 2p and 
_D 2 i/ > 9 (13. 2p . and repeatedly using the Cauchy-Schwarz inequality gives us 

> -C - C\D 2 u\ 2 + fi \--f\D 2 u\ 2 -C 



2 

+r l (5"\D 2 uDH\ 2 ( 1 \ - r]\D 2 u\) + 

+P' I A Q^|£> 2 m| 3 - C - C|L> 2 «|) - Crj\D 2 u\ 2 - CVP 2 m| 3 J . 

Here C is denotes various constants that are independent of e, but may depend on \DH(Du(xq))\, 
t](x ), Dt](x ), D 2 7](x ) and 

C := max| J D 2 if(Dw(p))| > (D^Du^o))!- 

Observe that C is a universal constant by assumption (13. 211 . 
3. Set 

A = A e := aM t 

where a > max{l/7,C"} is a positive constant. Notice that 

">2i2 , ,, ( 1,1 n2„,|2 



> -C - C\D'u\ 2 + n [ --f\D 2 u\ z -C 



+7 1 2 \D 2 u\P"\D 2 uDH\ 2 { 1 a - 1) + 

+P' ^ar]\D 2 u\ \^ m 6\D 2 u\ 2 -C- C\D 2 u\ \ - Cri\D 2 u\ 2 - C'r] 2 \D 2 u\ 3 ^ 



\D 2 u\ 2 + fi[^\D 2 u\ 2 -C 



+P' ^ar]\D 2 u\ (^ri~f9\D 2 u\ 



2 -' 2 C - C\D 2 u\ - C V \D 2 u\ 2 - C'r] 2 \D 2 u\ 3 } . (3.9) 



If for this choice of a 



ar]\D 2 u\ Qt770|D 2 «| 2 -C- C\D 2 u\^j - Crj\D 2 u\ 2 - C'r} 2 \D 2 u\ 3 < 0, 

then there is a bound on rj(x ) \D 2 u(x ) \ depending only on 9, \DH(Du(x ))\, and C. If the 
above inequality does not hold, then (13. 9p implies 

l2„.|2 1 ..(1,1 r>2„.|2 



> -C - C\D 2 u\ z + fi [ -i\D l u\ l -C). 
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For /i chosen large enough, this inequality implies a universal bound on t](xo)\D 2 u(xq)\. 

In all cases, we have bounded r](xo)\D 2 u(x )\ 2 from above independently of < e < 6 
(but perhaps dependent on 8, and C) and therefore, 

M 2 < irmxv < rj(x Q )\D 2 u(x )\ 2 + CM, + C < C(M e + 1). 
o 

Consequently, M e < Ci for some Ci (possibly) depending only on O', C and 6. □ 

4 Convergence 

In this final section, we complete the proof Theorem 11.11 We have already established the 
existence of a unique viscosity solution in Proposition 12. 1( see Remark 12.21 The only issue 
remaining is the regularity of solutions, which will be handled by our results in the previous 
section. Note however, that we established Proposition 13 . 1 1 under the extra hypotheses ( 13. 2ft . 
while Theorem 11.11 does make this assumption. The main idea to overcome this technicality 
is to employ the inf- convolution 

H\p):= inf lH(q) + ±-\p-q\ 2 ), (p, t) G W 1 x (0, oo). (4.1) 

We make use following lemma and omit the proof as it is elementary. 

Lemma 4.1. Assume that H : M. n — > R satisfies ( II. 3p . 

(i) H\0) < 0, fort> 0. 

(ii) H l H locally uniformly onM. n , as t — > + . 
(Hi) H l is convex and for each p, z el B and t > 

\z\ 2 

H\p + z)- 2H\p) + H\p - < — . 
In particular, < D 2 H f < 1/t, a.e. peM". 

Assume that H satisfies (jl.4p and D 2 H(p) > 6, for a.e. p G M n . 

(iv) For each t,R > 0, 

H\p + z)- 2H\p) + H\p -z)< (\D 2 H\ L ^ iBQ{R t)+lzl) ) \z\ 2 
for \p\ < R and z G M n , where 

Q(R,t) 

(v) For each p, z, G R n , t > 

H\p + z) 

Hence, D 2 H l > 9/(1 + t9), a.e. p G R n , t > 0. 



2 ( R + t max \DH(w) 

\w\<R 



2H\p) + H\p- z) > 



e\z 



2 



1 4- ±0 



1(3 



Proof, (of Theorem ll.il) 1. First assume that H satisfies (13.21) . so that we can apply Propo- 
sition 13.11 In view of the conclusion of this proposition, a routine application of the Arzela- 
Ascoli Theorem implies there is a sequence of e k — > such that 

u tk — > v uniformly in O 
u^^v inCt(O) 

as k — > oo. That is, u €k — > v uniformly in O and u tk — > v in C l (0') for each O' GG O, as 
& — > oo. It is clear from the above convergence that v G C{£(0) n C°'\0). We now claim: 
u is a viscosity solution of (jl.ip and therefore has to coincide with u by the uniqueness of 
viscosity solutions of (11. ip . 

2. Suppose that v — ip has a local maximum at Xo G O and that </? G C 2 (0). We must 
show 

max {5v(x ) - a(x ) ■ D 2 <p(x ) - f(x ), H(D<p(x ))} < 0. (4.2) 

By adding x >->■ f |x — x | 2 to and later sending p — > 0, we may assume that v — ip has 
a s£n'ct local maximum. Since u tk converges to v uniformly (for some sequence e k — > 0) as 
k — )• oo, there is a sequence of Xk such that 

Xk — > %o, as k — > oo 

— </) has a local maximum at x k 

As u £fe is a smooth solution of ( 11.51) . we have 

5n efc (x fc ) - a(x fc ) • £>V(xfc) + Pe k (H(D<p{x k ))) < f(x k ). 
Since /? £ > 0, we can send k — > oo to arrive at 

Sv(x ) - a(x ) ■ D 2 <p(x ) < f(x ). 

By Lemma [3.51 

< P £k (H(Dcp(x k ))) = !3 ek (H(Du' k (x k ))) < C, 
which necessarily implies that when k — > oo 

H(Dip(x )) < 0. 

3. Now suppose that v — ip has a local minimum at Xq G O and that ip G C 2 (0). We 
must show 

max{6v(x )-a(x )-I?rl>(x )-f(x ),H(DrJ;(x ))}>0. (4.3) 
Arguing as above, we discover there is a sequence e k — > as k — > oo, and such that 

%k — > x o, as k — 7- oo 

— ^ has a local minimum at x k 
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If 

H(Dt/j(x )) > 0, 

then (14. 3p holds. Suppose now that 

H(D4j(x )) < 0. 

Since u e is a smooth solution of (II .5j) . we have 

Su ek (x k ) - a(x fc ) • Dfy(z fc ) + ft t (F(^(i t ))) - /(z*) > 0. (4.4) 

By the convergence established in part 1 of this proof, H(Difj(xk)) = H(Du ek (xk)) < for 
all large enough k. Hence, 

lim p ek (H(Lh/>(x k ))) = 0. 
In this case, the above limit and (I4.4p imply 

max{5v(x ) - a(x ) ■ D 2 tf}(x ) - f(x ),H(D^(xo))} > 5v(x Q ) -a(x ) -D 2 i/;(xo) -f(x ) > 0. 

This argument verifies the claim that v = u and in particular, it proves Theorem 11.11 under 
the hypothesis that H satisfies (13. 2ft . 

4. Now assume that H satisfies (II. 3p and define for t, p, G (0,1) the function 

H t,p,0 (p) :=6\p\ 2 + p G M n , 

where 

H^(p) := / ^"(p - y)H t (y)dy, p G E n . 

Here if* is the inf-convolution of H given in (14. ip . t] p (x) = p~ n r](x/p), rj G C%°(M. n ), 
J Rn rj(x)dx = 1, and rj(x) = for \x\ > 1. Therefore, G C 00 (M ri ), and as p tends to 
0, converges to H l locally uniformly (consult Appendix C in [3] for more on mollifica- 
tion) . 

By parts (i) and (in) of Lemma [4.11 H t)P) satisfies (13.21) and (for p small enough) (II. 3p . 
Hence, 

J max {Lu - f , H*'* 6 \Du)} = 0, x G O 
\ u = 0, x G dO 

has a unique solution u t,p,e G C^. Moreover, an easy conclusion of the convergence assertion 
in parts 1 — 3 of the current proof gives the estimate 

lu^V-co') <C(0',a). 

for each t, p,9,a G (0, 1) and O' CC O. It follows that there is are sequences pk, 9k — > 
and u tk ' pk ' 6k converging in Cl oc (0) to some u G C^(0) fl C 0,1 (O), as A; — )■ oo. By part (ii) 
of Lemma 14.11 and the stability of viscosity solutions, u is the unique viscosity solution of 
(11. ip . This completes the proof of part (i) of Theorem 11.11 
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5. Now assume in addition that H satisfies (jl.4p . According to parts (iv) and (v) of 
Lemma |4TT| H t,p,d also satisfies (jl.4p . By assumption, D 2 H > 9' for some 9' G (0, 1). Lemma 

implies that D 2 H l ^ e > \9' , for t > small enough. The convergence assertion in parts 

1 — 3 of the current proof implies the estimate 



\ r ,, t -Pfi\ 



<C( O',0',iiiax|D 2 fr*'^(p)| ] 
V \v\< c J 



for t,p, 9 > and some universal constant C. Noting that part (iv) of Lemma [4.11 has that 
D 2 H t ' p ' e is a.e. locally bounded above, independently of all t, p, 9 positive and small, we are 
able to conclude as in part 4 of this proof. □ 

As a final remark, we mention equations with general boundary conditions 

max{Lu — f, H(Du)} = 0, xEO (4:5) 
u — g, x G dO 

can be handled similarly to PDE (11.11) . The method presented in this paper works with little 
alteration provided there exists a subsolution u G C x (0) of (14. 5 p such that u\qq = g and 

H(Du(x)) < 0, x GO. 
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